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Abstract 



In this paper we study the influence of the threshold effects due to massive 
degrees of freedom in the evolution with scale of gauge coupling constants. 
We first describe in detail the (standard) mass dependent renormalization 
prescription we use. This guides us to introduce and work with effective cou- 
plings, which are finite, process independent, and include complete thresh- 
old effects. We compute the evolution of the effective couplings in both, 
the Standard Model and its Minimal Super symmetric extension, from mz to 
the high energy region. We find that the effects from thresholds due to the 
standard massive gauge bosons are non-negligible, contrary to what is gener- 
ally assumed when using other, less-accurate, descriptions of the thresholds, 
as for example in the step-function approximation. Moreover, we find that 
thresholds are relevant when studying perturbative SUSY unification, chang- 
ing the conclusions reached when using the step-function approach. We find 
that threshold effects bring conflict between the known experimental data at 
mz, the naturalness upper bound on the masses of SUSY partners and the 
perturbative unification of couplings. 
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I. INTRODUCTION. 



In studying Grand Unification theories for the strong, weak and electromagnetig interac- 
tions, we are compelled to work with two different scales: the electroweak symmetry breaking 
scale (0(mz)), on one hand, and the scale of the symmetry breaking of the unifying group 
G, or unification scale (Mx ~ 10 ie GeV), on the other. In general, the unification schemes 
assume the existence of a "desert" between these scales (separated, at least, by 13 orders of 
magnitude), so that no new effects are present, beyond those coming from the unification 
group G. 

To test the validity of the unification picture, one studies the evolution of the three 
gauge couplings of the Standard Model, from low (mz) to high energies (Mx), seeking for 
their convergence towards a common value at and beyond the unification scale Mx- In the 
simplest scenario, we can neglect the effects of the heavy degrees of freedom introduced by 
G in the evolution of the running couplings, and work only with the content of matter of the 
low energy theory (light degrees of freedom) . This would be a first step in towards checking 
unification for the model. If, in this approximation, the couplings meet at a common value, 
g, at some point of the high energy region, Mx, we will have given the first step in the right 
direction. 

The following step would consist of including the effects of the heavy degrees of freedom in 
our scheme. At low energy, these degrees of freedom are decoupled from the theory, but they 
become operative as soon as we approach the scale Mx, changing smoothly the evolution of 
the couplings. To take into account their effects in a proper way, one can integrate out the 
heavy fields from the complete action of the full theory, S[G] Carrying this out, one 

obtains as unification condition: g~ 2 (rn? = g~ 2 (/i) + Aj(/i, Mj), where the functions Aj(/i, Mj) 

depend logarithmically on the masses of the heavy fields Mj rs O(Mx), the scale /i satisfies 
the condition m, <C \i <C Mj, with m; the masses of the light fields. With this construction, 
the behavior of gi(fi) is only dictated by the light degrees of freedom. 

The main tool to carry out this two-step program are the renormalization group equa- 
tions (RGE) and in particular the /3- functions for the coupling constants, defined as 

A = 77 — ■ in most of the Grand Unification analyses carried out, the /3-functions used 
a In /i 

are those computed using the modified minimal subtraction (MS) scheme ||. But, since 
there is no unique renormalization prescription to work with, and different prescriptions 
yield different explicit forms of the (5- functions, the question that naturally arises is to 
which extent the results on unification depend on the prescription chosen, that is, on the 
way one studies the evolution of the coupling constants with scale. 

The choice of renormalization prescription is not a matter of taste. It implies the way one 
deals with physical effects, such as threshold effects, coming from the presence of massive 
particles in the theory. The MS procedure is only one of the possible choices of a generic class 
of subtraction procedures, called mass independent subtraction procedures (MISP) [|J. For 
this class, the /3-functions depend only on the particle content at the energy scale at which 
one computes them, and not explicitly on the masses of the particles. On the other hand, 
there also is available another generic class of procedures, called mass dependent subtraction 
procedures (MDSP) 0, which take into account the dependence not only on the number, 
but also on the particle masses. While with the resulting MISP-/3-functions one is forced 
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to put in by hand the information about the mass spectrum of the theory when integrating 
the RGEs, with MDSPs one includes, in a natural way, all the information about threshold 
effects and possible decoupling of the massive particles in the evolution of the coupling 
constant ||. Because of this feature, the MDSPs are threfore more complete and precise 
than the MISPs. They have the disadvantage that the calculation of the /3-functions in 
MDSPs is far more complicated from a technical point of view. 

Working with the MS /3-functions, the standard procedure to take into account the 
contribution of a particle of a given mass m at same scale /i, consists in making use of 
a "step-function" 8(fi 2 — m 2 ): one for each massive particle J7|. This method constrains 
the contribution of a particle strictly to scales higher than its mass, being zero otherwise, 
and gives only the dominant logarithmic contribution of the mass. On the contrary, in the 
MDSP-/3-function, the contribution of each massive degree of freedom is controlled via a 
smooth function /(m//i), calculated in perturbation theory, that has the limits f{m/ji) — > 1 
when the ratio (m/fi) goes to zero, and f(m/fi) — > when (m/fi) goes to infinity. They 
give a non-zero contribution even for scales /x < m, and a more accurate description of the 
threshold crossing than the step-function approximation (see for example, Fig. fllCf)). This 
behavior is a quantum-mechanical effect, reflecting the Heisenberg uncertainty principle: we 
have a non-zero probability of producing a particle even for momenta below the mass scale 
of the particle; furthermore, the contribution of the degree of freedom to the /3-function 
spreads over a few orders of magnitude in momentum. 

On the different treatment of threshold crossing resides the main source for the discrep- 
ancies in the final results derived with one or another approach. The change in the derivative 
of the coupling constant when crossing the threshold will affect the value of the couplings 
even when we are away from the threshold region at several orders of magnitude above the 
particle mass scale. And this will have clear implications in the study of unification theories 
with a rich low energy mass spectrum, such as for the Minimal Supersymmetric extension 
of the Standard Model (MSSM). 

Previous work on the unification of the minimal susy model indicates that unification 
of gauge couplings is possible, and compatible with a susy spectrum below 1 TeV and a 
unification scale of O(10 16 GeV) These results were carried out by using the MS 

/3-functions, corrected with step-functions to treat the susy spectrum. Here, we will instead 
work with a different renormalization prescription, a mass dependent prescription, to include 
threshold effects. We will compare the results with those obtained by other authors. In this 
way, we will try to understand where and why are the differences between the conclusions 
reached with different approaches. 

With the derivation we use to get threshold effects, we aim to clarify the point that these 
effects are not higher order corrections with respect the order of perturbation theory one is 
working, but instead they have to be calculated at each order of perturbation theory. 

Before dealing with the susy model, we need, of course, to fix our prescription. We 
would have to choose one of the possible MDSPs. However, better than pick up a particular 
procedure and compute the /3-function, we will see that the use of an MDSP allows us to 
define effective couplings, analogous to the effective charges defined by Gell-Mann and Low, 
Stuckelberg and Peterman in QED fllQfl , in the sense that they are finite, universal, and 
include complete threshold effects, being by the way appropriate for our purposes. This will 
be the subject of Section 2. We will end that section with the definition of gauge effective 
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couplings associated with the gauge group SU (3) c x SU (2) l x U(1)y- Some problems related 
with the definition of the effective couplings for non-abelian theories, and also with broken 
symmetries, are also discussed in Section 2. 

In Section 3 we first study the evolution of the gauge couplings of the SM with minimal 
matter contents: three generations of fermions, and one doublet of scalars. In this case, the 
number of massive degrees of freedom in the range (mz, Mx) is small, and we do not expect 
the unification of couplings to take place. In spite of this, the model becomes useful to point 
out the importance of threshold effects associated with massive gauge bosons, W and Z°, 
otherwise neglected in the step-function approximation. 

We should keep in mind these contributions when studying the minimal susy extension 
of the SM, where more degrees of freedom (susy partners) are present beyond the scale mz- 
We will see in Section 4 that gauge boson threshold effects become relevant in attaining or 
not the "unification" of the gauge couplings. We first examine the model in the simplest 
approach (one-loop without unification gauge group thresholds) ; after that we will consider 
corrections due to heavy threshold effects in the standard way described in the beginning of 
this Introduction, and also will discuss the effects due to two-loop corrections. 

In Section 5 we offer our conclusions, and discuss other approaches that have emerged in 
the recent literature [Q-Q|. The works cited on Ref. [I3| remark the importance of complete 
light (susy) threshold effects for the determination of the MS values of the gauge couplings 
at the electroweak scale, as they are extracted from experimental quantities, and the later 
implication of these initial values for unification. Susy thresholds distinguish between the 
quantities extracted using the SM, from the same quantities obtained on the MSSM. The 
main point of these papers is that not only the leading (logarithmic) threshold contributions 
are important, but also the non-leading, non logarithmic terms. In some sense, that is 
also our main statement. But we remark that here we study the impact of thresholds in 
the evolution with scale. Since we aim at comparing different evolutions, as given by two 
different approaches for the thresholds, we will assume the same initial conditions for the 
couplings in both schemes^. 



II. EFFECTIVE COUPLINGS. 

Our purpose in this Section is to define effective couplings, following the simplest possible 
arguments from the general theory of renormalization. This is an easy task when the coupling 
is the coupling constant for an abelian gauge theory, as in QED. We will use what we learn 
in this case as a guide to get at the end the expression for non-abelian effective couplings. 

For a general theory, the renormalized coupling, g r , is related to the bare coupling, go, 

1 /2 1 

by g r = Zg , where Z is a product of renormalization constants: Z = Z 3 (Z 2 Z^ ) a . Here, 
Z$ is the wave function renormalization constant (WFRC) of the gauge boson which de- 
fines the interaction; Z 2 is the product of WFRC of the external legs; and Z\ is the proper 
vertex renormalization constant. The subindex "a" refers to the different classes of vertices 



Further work in the determination of the initial values for the effective couplings at the elec- 
troweak scale in our scheme, is now in progress. 
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involving the same coupling constant, which come from the interaction of the gauge boson 
with itself or with the other particles present in the model: fermions, scalars, .... When 
using a MISP, the Slavnov-Taylor identities (Ward-Takahashi for abelian theories) guar- 
antee that (Z 2 Z{ 1 ) a is independent of vertex we choose to define the coupling, and thus 
we obtain in this case a universal renormalized coupling, albeit without including thresh- 
old effects. With a MDSP, we include the dependence on the masses, in particular on the 
masses of the external legs, and we could, in principle, distinguish among the renormalized 
couplings derived from different processes: we would get a mass-dependent coupling, but 
no a universal coupling [13]. This would not be inconsistent with the constraints imposed 



by the Slavnov-Taylor identities. For example, one can define the renormalization constants 
associated with a specific process, and calculate the /3-function, and after make use of the 
Slavnov-Taylor identities to derive the remaining renormalization constants for the theory, 
without modifying the previously defined /3-functionP]. But the renormalized coupling would 
not be universal (independent of the process). Therefore, in order to get a universal renor- 
malized coupling, we first have to make sure that the MDSP respects the Slavnov-Taylor 
identities, using them to select the finite contributions (which contain the threshold effects) 
to the Z{. This will guarantee the gauge invariance of our effective couplings. 

For abelian gauge theories, such as QED, one can extend the Ward-Takahashi identity, 
Z 2 = Zi, to the finite terms of the Zj, and therefore the relation e r = Z^ 2 g is valid for 
mass independent and mass dependent procedures. Moreover, we can see that, if we choose 
a suitable mass dependent procedure, the renormalized charge is equivalent to the effective 
charge of Gell-Mann and Low, Stuckelberg and Peterman, defined by: 

Here, LTq (g 2 ) is the transverse component of the bare vacuum polarization tensor of the 
gauge boson. The effective charge is (a) independent of the renormalization scheme, since 
it is defined via bare functions, (b) universal and gauge independent and, (c) finite. 

The last point is easily demostrated when one replaces the bare coupling by the renor- 
malized coupling in (|2.1| ), since Z3 must verify that the product Z 3 (l + I1q (g 2 )) be a finite 
function of both the scale q and the subtraction point fi. We therefore have: 

;i+iw,/i 2 ))' 

where the explicit form of rpj(g 2 ,/i 2 ) will depend on the chosen subtraction procedure. 

In general, the two couplings e\^ and e 2 are not equivalent. For example, the renormal- 
ized charge using MS does not include threshold effects due to massive fermions, whereas 
they are present in e\^ through the functions n^g 2 ,^ 2 ). On the contrary, if we adopt 
a mass dependent procedure with the normalization condition 1 + ILj(/i 2 ,/i 2 ) = 1, then 

e lfM 2 ) = 4(^ 2 )i and 

e eff e r 



e efM ) = n i TT7Y„2 . .2\\ > ( 2 - 2 ) 



2 This was explicitly shown for QCD in Ref. 
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Therefore, both constants depend on the scale in the same way, and they include the thresh- 
old effects. 

In the case of non-abelian theories, to get a finite expression for the effective couplings it 
is necessary something beyond the correction to the vacuum polarization. We must supply 
another function, T (q 2 ), to the definition fl2.1|), 

9r(v) 

yeffKH ' (i+n V I M 2 ))(i + 2r ( ? 2 )) z 3 (i + uUq 2 ,^))(z 2 z^)i(i + 2r ( q 2 )) ' 

(2.3) 

which comes from the corrections due to vertex and external legs, and which must verify 
that the product (Z 2 Z{ 1 ) a (l + ^o(q 2 )) be a finite function. The function T (q 2 ) consists of a 
divergent term, Tf ,v , and a finite contribution which includes a dependence on the masses. 
The Slavnov-Taylor identities guarantee that the divergent term of (Z 2 Zf 1 ) a is the same 
for all possible vertices "a", and is given by fllEfl : 

( W)f = 1 + (J^(G)^±i , (2.4) 

where C 2 {G) is the quadratic Casimir of the non-abelian group G, £ is the gauge parameter, 
and n is the dimension of the space-time. Thus, 1 — Tf ,v = (Z 2 Z{ )f v . 

But, it is not straightforward to calculate the finite contributions to r (g 2 ), i.e., to 
define (Z 2 Z{ 1 ) a , with a mass-dependent procedure. The problem is to make sure that the 
resulting To(q 2 ) is independent of the vertex "a" chosen to perform the calculation. For 
example, if we choose the vertex with fermions on the external legs, we can expect that 
(Z 2 Z{ 1 ) f depends on the fermion masses, the scalar masses, and the gauge boson masses (if 
the symmetry is broken). The same kind of masses may appear if we consider the scalar- 
gauge boson vertex, or the trilinear boson vertex. However, if we choose the ghost-gauge 
boson vertex, we only have running in the loop, ghosts, gauge bosons or Goldstone bosons 
(if the symmetry is broken): either massless particles, for unbroken symmetry, or particles 
with masses proportional to the gauge boson mass, for broken symmetry. Now, if we impose 
the Slavnov-Taylor identities, and thus (Z 2 Z{ 1 )f = (Z 2 Zi l ) g h os t = . . ., it is clear that the 
universal correction To(q 2 ) can not depend on the fermion or scalar masses, but only on 
the gauge boson masses. 

In this sense, the universal threshold effects associated with the vertex are related to 
processes leading to gauge bosons production. Therefore, if we have an unbroken symmetry, 
there are no threshold effects from the vertex in the effective couplings, because we can 
always produce a massless particle. The corrections due to the presence in the vertex of 
other massive particles, such as fermions or scalars, will affect other parameters of the 
theory, but not the effective couplings. Moreover, these process depending corrections are 
finite, as it is implied by the Slavnov-Taylor identities. 

So far, we have discussed about the general expressions of the effective couplings (for a 
general non abelian theory) given in terms of the functions ITq (q 2 ) and T (q 2 ). The transverse 
bare vacuum polarization is easily calculated from the appropriate Feynman diagrams. We 
know the divergent part and the masses presents in the finite term of r (q' 2 ), but the explicit 
form of r (g 2 ) will be closely related to the specific gauge theory we treat. 



„2 („2\ 
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In particular, we are interested in studying the evolution of the gauge couplings of the 
Standard Model, g 3 , g and g' of SU(3) C x SU{2) L x U(l)y- We have no problems with the 
definition of g^q 2 )-, since QCD is an exact non abelian theory, so that: 



where, 



W) = -gh(S + - In \ + Constant , (2.6) 



) 



and IlJ 9 (g 2 ) is the transverse bare vacuum polarization of the gluon, given by the diagrams 
of Fig. 

To define the other two effective couplings, g and g', we have to take into account that 
SU(2) L x U(l)y is a broken symmetry at low energy (electroweak scale), just where we begin 
to run the couplings. In the broken phase, we have the three gauge bosons (the eigenstates 
of the mass matrix) W^, Z° and the photon, A. We can define the effective coupling g 2 (q 2 ) 
by the interaction of the W ± ; however, since in this phase the gauge boson B, associated 
with the gauge symmetry U(1)y, is a mixed state of the neutral gauge bosons, it is better 
to work with the electromagnetic coupling, e 2 (q 2 ), given by the interaction of the photon, 
and define g' through the equation: 

111 . , 

—2=---- 2 ' 7 

gU e Z gZ 

In the Standard Model, e 2 (q 2 ) is not a pure abelian coupling, so that to define g(q 2 ) 
and also e 2 (q 2 ) we need the correction T(q 2 , m 2 ), where now "m" may be mz or mw This 
vertex correction is the same for both couplings, because the non abelian character of e 2 is 
closely related to g 2 (e~ 2 = g~ 2 + g'~ 2 ). To obtain its explicit form we follow the argument 



given by Kennedy and Lynn [17]], and that we partially reproduced in the following. 

These authors make use of the relation between the non-abelian vertex correction, T(q 2 ), 
and the longitudinal term of the mixed vacuum polarization tensor for the neutral bosons, 
U ZA (q 2 ) [18]. The latter contribution gives rise to a non-diagonal mass matrix for the 



(Z°, A) system. Therefore, one needs to redefine the fields Z° and A, which are the correct 
mass eigenstates at tree-level, in order to eliminate the non-diagonal term, and get the 
correct eigenstates (and a massless photon!) at one-loop order. 

This can be carried out by first redefining the coupling g, including the universal vertex 
correction T with g = g(l — T), and then defining new fields Z and A making use of this 
coupling. Now, in this basis the non-diagonal term of the mass matrix becomes Hz a + 
m 2 z gg'T j \g 2 + g' 2 ). Thus, if we choose Y to satisfy the condition 

ni4 + mi-^r = o, (2.8) 

gl + gl2 

we get the desired results that: (a) the photon remains massless (at least at one-loop), and 
(b) we determine the explicit form of T(q 2 ). Obviously, condition (|2.8|) is an identity for the 
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divergent terms of Hza and T. What we get with (|2.8|) is the finite and mass-dependent 
term of T. 

Now, we have all the ingredients involved in the definitions of g 2 {q 2 ) and e 2 (q 2 ), which 
are given by: 

5V) = i + n T w (? 2 ) + 2r( 9 r e V) = i + uUf) + 2# rv) ' (2 ' 9) 

where n^^, U^ A are the transverse bare vacuum polarization tensor of the W ± and the 
photon respectively. 

In Appendix A we give the expressions of the functions 11^^, 11^ and T(q 2 ), as well 
as ITJ . All these functions, including r 3 (g 2 ), depend on the gauge parameter £. This 
dependence cancels exactly for the divergent part of the combinations LT^ + Tj, but the 
same does not happen for the finite term due to the presence of degrees of freedom whose 
masses are proportional to £ (i.e., gauge bosons, ghosts, and except in lTj g and r 3 (g 2 ), also 
Goldstone bosons). This leads to effective couplings depending on the gauge parameter |fTo| . 

The gauge parameter also changes with scale, as given by its RGE 

-e^— £, (2.10) 



d In q 2 d In q 2 

and it will depend on the associated^ coupling g 2 (q 2 ) through the dependence of Z 3 . In 
order to get the evolution of the effective couplings we have to solve a system of coupled 
differential equations. The change in gf{q 2 ) produced by the change in £(g 2 ) will be small 
(order two-loop); and if we calculate the effective coupling at one-loop order, we can neglect 
it, and maintain £(g 2 ) at its initial value. But if we calculate gf{q 2 ) at two-loop order we 
need, at least, £(<? 2 ) at one-loop order, and so on. On the other hand, the theory by itself 
indicates the most suitable value of £ to work with, without approximations: the fixed point 
of the differential equation at £ = 0. Thus, we choose to work in the Landau gauge. 

Up to now, we have defined the effective couplings in the Standard Model, which are 
finite, universal and include threshold effects. For the SU{2)l x U(1)y couplings, we made 
partially use of the arguments given by Kennedy and Lynn to derive the function T(q 2 ). 
Nevertheless, our effective couplings do not coincide exactly with their definitions. In par- 
ticular, we differ in the definition adopted for g(q 2 ). Here, we have chosen to relate this 
coupling directly to the W ± propagator, in the same way that the coupling e(q 2 ) is related 
to the photon propagator. They choose, instead, the mixed propagator of the photon and 
the W 3 boson. Their convention follows from imposing that the tree level relation among 
e 2 , g 2 and the sine of the mixing angle, sin 2 9w = s 2 ^, had to be maintained at the level of 
the effective parameters, i.e, 

4rWALtf) =1 (2 . n) 

e z {q z ) 



3 We will have four different gauge parameters, associated with the bosons W^, A and Z° and 
gluons respectively. 



S 



where s 2 v {q 2 ) is derived from neutral current processes. 

In our case, to get a definition of s^(g 2 ) consistent with the neutral current amplitude, 
we have to allow for the relation ( |2.11| ) to receive radiative corrections (T^j, so that, 



e 2 (q 2 



0(h). (2.12) 



The coupling g^iiq 2 ) is related to a conserved current, and hence it only receives contribu- 
tions from loops of charged and degenerate particles. There is, for example, no contributions 
due to the Higgs. On the other hand, all the doublets under SU (2) L contribute to the cou- 
pling g 2 (q 2 ), and this includes the Higgs. The coupling g 2 (q 2 ) is related with a broken 
symmetry, and this fact is not reflected in g 2 KL (q 2 ). We consider that the definition of g 2 (q 2 ) 
is more appropriate for our purposes, whereas gxiiq 2 ) is more appropriate to study the 
neutral current and related processes. 

In the following sections, we will proceed to study the evolution with scale of the effective 
couplings from mz to the high energy region, where we want to check for the validity of the 
unification scenario. With this in mind, we see that the above definitions of the effective 
couplings will be useful once we eliminate the bare couplings in favor of the couplings given at 
the scale mz- Furthermore, bearing in mind that the combinations IT? + 2rj are proportional 
to g 2 , we redefine 11^ + 21^ = ^f(II^ + 21^), and obtain: 

+ (Tlliq 2 ) + 2T % {q 2 ) ~ n£(m|) - 2T t (m 2 z )) . (2.13) 



gf{q 2 ) gK m i 

As an example of unifying group we will take the minimal choice, i.e., G = SU(5). 



III. THRESHOLD EFFECTS IN THE STANDARD MODEL. 

In this section we will study the Standard Model with minimal matter contents, i.e., three 
generations of fermions and one scalar doublet. From the point of view of unification, we do 
not expect that the threshold effects included in the definition of the effective couplings can 
change the negative results obtained with other renormalization procedures M . This model 



is simply intended to see how thresholds affect the evolution of gauge couplings [2(| . 

Since for scales fi > mz we can regard all the fermions as being massless, except for the 
top quark, the relevant masses for the problem are mz, mw, mt (top) and nih (Higgs). The 
values of mz and m w are determined by experiments to be at pi| : m z = Ol.lSTiO.OOTGeV^, 
mw = 80.2 ± 0.3 GeV. But m t and mh remain as free parameters, with lower experimental 
bounds^ m t > 131 GeV, mh > 60GeV [^]. On the other hand, imposing as condition 
the validity of perturbation theory in the range of scales considered, one gets the upper 



The recent data from the CDF Collaboration at Fermilab [22] indicate the existence of the top 
quark with mt 



174 ± lOljJ GeV . However, this data is still to be confirmed, and therefore we 



would maintain the lower bound given in Ref. [21]. In any way, to fix or let the top mass free makes 
only some very tiny differences in the numerical results, and it does not influence the conclusions 
reached in this paper. 
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theoretical limit: m t , rrih < 200 GeV |[24|| . We will use in this section these latter values for 
m t and m^. 

The initial values of the couplings, g 2 (mz), are identified with the experimental values 
measured at the scale m z : 

a 3 (m z ) = 0.117 ±0.005, [21] (3.1) 
a~\m z ) = 127.9 ±0.3 m (3.2) 



sin 2 9w(mz) = 0.2319 ± 0.0005 pl|. (3.3) 



In what follows, we will work with the "constants" a,{ = gf/An. The initial value of a^ 1 is 
fixed by the relationship: (mz) = a~ x {mz) sin 2 Qw( m z)- And, as mentioned above, the 
evolution of ct{ X is calculated by the relation a^ 1 = 3(a~ 1 — a 2 ^ 1 )/5, where the normalization 
factor is due to the embedding of U(l)y in 577(5). 

In Fig. (^) we have plotted the effective couplings (a" 1 \ e f) in the high energy region. We 
have also plotted the couplings calculated in the step-function approximation (c^" 1 \ e ), and 
using minimal subtraction (a" 1 \ms)- The latter procedure is equivalent to the step-function 
approximation when all masses are below or at the scale mz- We see in the plot that the 
approximation we use for aj -1 and a 3 l makes only very little difference. Nevertheless, we can 
see that, due to threshold effects produced by the top, a 3 1 \ e f is a little larger than a^ 1 \ e , 
and this, in term, a little larger than a 3 l \ms- The top quark decouples at low energies 
{m t <C 10 12 GeV), but its effect propagates to the high energy region by the RGE, being the 
decoupling in a 3 l | e / smoother than in the ^-approximation. 

For the coupling a^ 1 , one would naively have expected that the threshold contributions 
due to m t and contribute to increase its slope. But these are not enough, because in this 
case the dominant threshold effects are those due to the massive gauge bosons, so that the 
slope of a^ 1 \ e f decreases, as can be seen from the plot. This effect is not present in a^ 1 \$ 
(and a^ 1 \ms) because we begin to integrate precisely at the scale of m z . This can be seen 
clearly if we compare the derivatives, given by: 

(47r) l0 = f ^ " m|) ~ T2 9{ ^ - ml) " - m2) - T2 ' (3 ' 4) 

(47T ) r = y( s w/s(0) a w) + Cwf g (a z , a w )) + 3f r (a w , a w ) 

1 37 

--f s (0,a h )-f f (0,a f )--, (3.5) 

where the functions fb{cii,aj) are defined in Appendix A, and smooth out the "step" of 
the Heaviside function, as demanded by the Heisenberg uncertainty principle (see Fig. [|). 
The derivative of a^ 1 \e changes from negative to positive for q 2 = m%, while in the case 
of a^ 1 \ e f threshold effects postpone this change of sign |TJJ]. Therefore, ce^ 1 \ e f decreases 



when we begin the integration, and this makes the differences with ct^ 1 \e sizable only in the 
high energy region. 

The differences between step-function and effective couplings are in that for the former, 
thresholds are considered as an "instantaneous" effect; we "switch on" and "switch off" them 
at a specific point of the energy scale, i. e., at q 2 = m 2 . But indeed, thresholds are spread 
over a certain range in the momentum scale. The threshold crossing is not "instantaneous". 
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As shown in Appendix A, only when we begin to integrate at a scale q$ ^> m 2 (massless 
approximation), or we end the integration at q 2 <C m 2 (complete decoupling), we get the 
same contribution for the massive degree of freedom in both a^ 1 \g and a^ 1 \ e ff. In the 
general case, when q 2 < m 2 < q 2 , even if we can apply the limit m 2 <C q 2 , and then suppress 
all the terms 0{m,i/q) in the threshold function contribution, we will get that, 

a i ~ In — , 

while, 

2 

a" 1 \ e ff hi-^- - lnCj , 

where Cj is a "constant" of 0(10)0. This constant contribution is not more than a reminder 
of the fact that, although if the particle is decoupled at q , in going to q we have to cross 
the region of scales O(mj), where the degree of freedom is neither decoupled nor coupled^. 
Moreover, the transition between the decoupled-coupled regimes, takes at least more than 
one order of magnitude in the energy scale. Therefore, when beginning the integration at 
the electroweak scale, mz, it is not a good approximation to consider the massive gauge 
bosons as completely coupled. 

We can now compare the corrections due to threshold effects with those coming from 
improving the order of perturbation theory in the mass-independent approach. In Fig. 
(|) we have plotted a^ 1 \ e f at one-loop order, and a" 1 \ms & t two-loop order. We also 
include a" 1 \ms & t one-loop order as a reference. We can now see that two-loop effects are 
qualitatively similar to threshold effects. Both of them raise a 3 1 and decrease a^ 1 , although 
for a? 1 the thresholds from the gauge bosons dominate over two-loop effects. We can expect 
that, when we take into account both corrections, calculating effective couplings at two loop 
order, the differences relative to the values of a^ 1 \ms become more accentuated. 

Instead of doing the exact calculation, we have decided to study this case in an approx- 
imate way. We use it only as an indication of which kind of behavior we can expect at 
2-loop order with thresholds. In Appendix B we describe in detail the argument followed to 
obtain approximated two-loop effective couplings, based on the expressions for the RGEs at 
2-loops, and general properties of the threshold functions. When we study their evolution, 
we find what we expected: a^ 1 decreases and a J 1 increases slightly, with respect to the one 
loop a^ 1 \ e f values. Of course, these effects are not enough to "unify" the three gauge cou- 
plings, but it is useful to keep them in mind when studying theories with a larger presence 
of massive degrees of freedom. 

In spite of the no unification of the couplings, the Standard Model with minimal matter 
contents has been useful to study the threshold effects introduced in the effective couplings. 



Constant" means that this term can be approximated by a constant. See Appendix A for 
details. 

6 The word "coupled" is used here as opposite to "decoupled", i. e., maximum contribution to the 
(3 functions. 
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In particular, we have seen that qualitatively, the thresholds associated with the gauge 
bosons are not negligible, contrary to what is assumed by the step-function approximation 
(since we begin to integrate at the scale q 2 = m z ). Quantitatively, moreover, this correction 
is larger than the two-loop correction. 

It is clear that thresholds are additive and they will become more relevant the larger the 
number of degrees of freedom in the range [mz,Mx]- This is the case of the Minimal Su- 
persymmetric extension of the Standard Model (MSSM), that we will study in the following 
section. 



IV. UNIFICATION IN THE MSSM. 

If supersymmetry is a symmetry of Nature, every particle of spin "j" of the Standard 
Model must come with its supersymmetric (susy) partner, of spin "j ± 1/2". Thus, we 
have gauginos, sleptons, squarks and higgsinos in addition to the content of matter in the 
Standard Model. We also need, at least, two Higgs doublets (with their higgsinos) in order 
(a) to cancel the anomaly due to the higgsinos, and (b) to give mass to the two components 
of the quark doublet. Since the susy particles have not been detected, susy must be broken 
at low energies, and therefore we have a larger number of arbitrary parameters of masses to 
consider beyond the scale m z . 

One can try to constraint the susy spectrum in the context of unification theories, making 
use not only of limits on proton decay p6|j2~T|| , but also of cosmological arguments on the 
relic abundance of the lightest susy particle (LSP) p8| . In general, the susy spectrum 
thus obtained lies below 1 TeV, and thus it is therefore compatible with the stability of 
the hierarchy of scales in the model under radiative corrections. In this sense, the limit of 
"1 TeV" is usually denoted as the "naturalness" bound. 

Most of the predictions have been obtained making use of the step-function approxima- 
tion to deal with susy thresholds. But the more accurate description of threshold effects, 
such as given in the effective couplings, could change the conclusions reached as we will show 



in this section [ 29| . 



Of the susy spectrum, we will assume that susy is broken at the Planck scale by a 
"hidden" sector, so that the large number of susy masses can be determined at the weak scale 
in terms of a small number of parameters given at the unification scale [[30]. Accordingly, 
we take the following parametrization, neglecting the (higher order) mixing between winos 
and higgsinos, and the s-tops left and right for simplicity: 

(a) For winos (w), and gluinos (g), 

= mi/2 , rn- g = 3m 1/2 , 

(b) for squarks (g), and sleptons left {II) and right {Ir), 



ni 



m l/2( C i + ft)) , C<j = 7 , Cj = 0.5 , Cf = 0.15 , 



where £o = (^o/ m i/2) 2 , wio being the common mass for the scalars, and rrii/2 being the 

gaugino mass at the unification scale; and 

(c) for higgsinos, m~ h will be taken as arbitrary. 
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Furthermore, there are five massive fields coming from the two Higgs doublets: two 
charged Higgses (m + ), two neutral Higgses {m h , m H ), and one pseudoscalar (m a ), whose 
masses must satisfy the relations: m 2 + = m 2 +my V , m 2 H = m 2 a + m 2 z — m\. The lighter neutral 
Higgs is identified with the standard Higgs, and for the remaining we take m 2 + ~ m 2 m m 2 H . 
At the end, we have as arbitrary mass parameters the following: mt, irih, mi/2, mo, m h an d 
rriH- The value of mi/2 is bounded from below by the experimental searches for charginos 
7L \, irii/2 > 45 GeV. We take the remaining susy parameters at least of order mz- 

The expressions of the effective couplings given in Section |HI|, with the functions Uf i + 2T i 
calculated in Appendix A, do not include the contribution of the susy degrees of freedom. 
However, these are easily derived, taking into account that (a) they do not appear in r i; but 
only in n^, and (b) they are contributions from fermions or scalars (Eqs. (|A11|) and (|A10|) ). 

As a first attempt to look for the unification of the effective couplings, we simply study 
their evolution with scale, without making any reference to the unifying group. In Fig. (f|.a) 
and (§[b) we have represented a" 1 | e / for two different values of the susy parameters (the 
lower and upper bound). In order to make a meaningful comparison, we have also included 
a" 1 \e calculated with the same susy parameters. Taking into account the experimental 
errors in the coupling constants, we see that the a~ l \g unify (they cut at one point, Mx, of 
the energy scale), for both values of the susy parameters. Due mainly to the experimental 
error in a 3 (m^), at this level of approximation one can not extract much more information 
about the susy spectrum with the a" 1 \g. 

But the situation is not as unambiguous for the \ e f. to get unification, effective 
couplings prefer higher values for the susy parameters, around the naturalness bound of 
1 TeV . Now, the threshold effects due to the gauge bosons play an important role in reaching 
or not unification. While the contribution of susy masses (fermions and scalars) tends to 
increase the slope of 1 \ e f and 1 \ e f, for the coupling a^ 1 \ e f this effect is compensated 
by the decrease of the slope due to the gauge thresholds. Because of this, the differences 
between a? 1 \ e f and a^ 1 \e are almost negligible, contrary to those obtained for the other 
two couplings, making a^ 1 \ e f move away from the "crossing" point with the other two 
effective couplings. Indeed, if the bosons W ± and Z° were massless, the slope of a^ 1 \ e f 
would increase sufficiently so as to allow the unification of the three couplings even for a 
susy spectrum of 0(mz) (see Fig. @). Qualitatively, we will achieve the same conclusions 
as we would if using the step-function approximation: the perturbative unification of the 
couplings is compatible with a susy spectrum below 1 TeV, and a unification scale Mx high 
enough to fulfill the experimental lower limit on proton decay. However, when we take into 
account the thresholds, including those of W and Z°, this general conclusion is slightly 
modified, and we get more information about the susy spectrum: it can not be of 0(mz)- 
We see that the electroweak breaking, being a "low energy" process (order mz), has effects 
on processes which take place at "high energy" (order of Mx), such as the perturbative 
unification of the couplings. 

Simply plotting the evolution of the couplings is not the best way for making predic- 
tions, and even less so with the big uncertainty we have in the experimental determination 
of as(rriz). Indeed, of the three gauge couplings the value of a^ijnz) is the worst known 
experimentally, in the sense that there is not a general agreement between different mea- 
surements |3T|-p3|. In general, the values obtained from experiments at low energies (below 
mz) are lower than the values inferred from the decay of the Z° into hadronic states (LEP 
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data) (see Table |). Because of this discrepancy it is better not to take any value of a^(mz) 
as the initial data and, instead, to derive it from the unification condition: 



a^(M x ) = a 2 \M x ) = a^(M x ) = c^ 1 . (4.1) 

In this way, susy masses can be bounded by demanding that 0.108 < 03 (mz) < 0.125, which 
are the experimental lower and upper values quoted in Table |. 

Solving the system of equations derived from ( |4.1|) , one gets the general expressions: 



I I 771 ■ 

<%\mz) \e = -(l^\m z ) - ^-\m z )) + — £ */(— ) , (4-2) 

7 OD7T j m z 

In— |» =^(3«^)-8«2 1 K))-iE</(-) 1 (4-3) 
mz 14 84 *-r mz 

where the sums run over all the masses (susy and non-susy). When we work within the 
step-function approximation, the functions "/" are simply given by ln(m 2 /m|); thus , only 
masses mj greater than mz contribute. Within the effective couplings approach, and making 
use of the approximation for the threshold function given in Appendix A (Eq. (??)), we get 

f{rm m z ) = In z f 1 , 4.4 
m| 

where the value of the q mostly depends on the kind of massive particle running in the loop 
(see Table 11). As we have seen in Sect. 3, and also in Appendix A, this constant is intended 
to smooth the threshold-crossing, and to control, not only beyond but also below, how far 
it is the threshold from the electroweak scale. 

Therefore, we see that, independently of whether we use effective couplings or the step- 
function approximation, the behavior with susy masses of (mz) and lnMj will be the 
same: the higher the susy masses, the higher (mz) and the lower is In M x . But now, with 
the effective couplings we have also a non zero contribution of the massive gauge bosons, 
which are explicitly given by: 



1 





-1.44 



L 2 i m 2 z + 2m w 2 m 2 z + 5m w 2 m| + 2.5m 
39s w In — \ — 7, 52 In 7, 36 In 5 — 

V m z + om^y m z m z 

o , m% + 2mw 2 , mi + 5mw 2 , m% + 2.5mw 2 1 
26gLln w , +651n z „ w +451n z - 9 — — ^+0.88. 

m z + 5m^/ m z m z J 

Thus, these terms decrease the value of a^ l (mz) and increase the value of InMx, with 
respect to those obtained in the step-function approximation. 

This behavior is clearly seen when we compare Fig. (0) (effective couplings) and Fig. @ 
(step-function), where we have plotted a^ l (mz) versus logio(mi/2), for different values of the 
remaining mass parameters. The allowed region in the map ctj 1 (mz) — logiom 1 / 2 is delimited 
by the experimental bounds on a 3 _1 (m^) and m\ii given above, and also by the theoretical 
limit M x > 10 16 GeV, which is a safe bound for proton decay. Since higgsinos and heavy 
Higgses contribute with the same sign, we have simply taken m~ h = m#. Moreover, a^ l (mz) 
and InMx depend only very slightly on the parameter £0 = (^o/' m i/2) 2 , so that we have 
taken £0 — 1 as a typical value for the plot. As we increase the mass of the higgsinos, we 
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decrease (increase) the value of a^{mz) (cnj i m z)) compatible with unification. However, 
since the step-function procedure does not distinguish masses lower than mz with any 
higgsino mass below mz and gauginos masses of 0(mz), we get the limit a^{mz) < 0.121, 
perfectly compatible with the experimental band for this coupling. As we have previously 
seen, no constraint on the susy spectrum is obtained at this order. On the other hand, 
with the effective couplings we would get the upper value of a 3 (mz) in the limit m~ h <C m Zl 
but, now, due to the contribution of the gauge bosons, this limit has no physical interest 
{oiz{mz) < 0.140). The point is that also for higgsinos masses of 0(mz), and almost for a 
susy spectrum of 0(1 TeV), the value we get for a^{mz) is too high. Turning the argument 
upside down, if we demand a^imz) < 0.125 and Mx > 10 16 GeV, we get a lower bound on 
= rriH, and an upper bound on mi/2- These bounds depend on m t and (the lower m t 
and nih, the higher m\/2 and m^), and also on £ , although these bounds remain practically 
unchanged for £ > 10 4 . For this value of £o we get: 

m 1/2 <7TeV (m h = 60 GeV , m t = 9lGeV), 
m~ h > 370GeV (m h = m t = 200 GeV , m~ h = m H ) . (4.7) 

By giving these bounds and the above bounds on a 3 (m^), we have taken into account the 
experimental errors in sin ^(m^) and a~ x {mz\ For the central values of these quantities 
we would obtain: rrii/2 < 3.8 TeV, m~ h > ITeV. Any kind of prediction from unification 
is very sensitive to small variations in the input parameters, and even small errors in the 
experimental data induce serious uncertainties^. 

With this calculation we have recovered in a more detailed way the general results we 
obtained before. At the one-loop order, we can get the perturbative unification of the 
couplings, compatible with the naturalness bound for the susy masses and proton decay- 
lifetime limits, independently of whether we take or not threshold effects into consideration. 
But the inclusion of more accurate thresholds than those described by the step-function, 
drives the susy spectrum close to the ITeV bound and the value of a 3 (m z ) into the range 
of LEP data. 

Threshold effects as considered here are not higher order corrections, but a correction to 
the simplest step-function approximation. As a first approach, just at one-loop order, this 
latter procedure gives a good indication for unification with a viable susy spectrum. Nev- 
ertheless, when one improves the treatment of thresholds, working with effective couplings, 
although the conclusions are not drastically changed, the result at one-loop does not leave 
much room for more corrections, for example 2-loop order corrections. 

In fact, it is not necessary to carry out a complete calculation of thresholds at two-loop 
to get the general behavior. It is enough to take into account that two loop effects tend 
to lower the value of a^ l {m z ) by about 10 %, for both the step- function and effective 
couplings^. In the first case, this effect drives the upper limit obtained at one-loop order 



7 Note that the values quoted in [29| are different, due to the different value of sin 2 9w{ m z) 
considered in that work. 

8 For the effective couplings, this can be easily shown by resorting to the kind of approximation 
given in Appendix B. 
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towards the experimental upper limit (0:3(777.2) < 0.125), while the naturalness bound on 
susy masses gives now the lower bound a 3 (m z ) > 0.121. With the effective couplings, a 
10% decrease pushes the values of a^ l {m z ) obtained with m~ h = m H < ITeV away from 
the experimental band. We would need masses of order lOTeV to get a^imz) < 0.125. For 
the two-loop calculation, the dependence on the parameter £0 is not so mild, and the value 
of a^ l {mz) increases with it. In any case, to have m~ h = mu in the range of TeV, we would 
need £0 > 10 10 , i.e., squarks and sleptons of order 1 PeV . Therefore, we can not conclude, at 
two-loop order, that the simple scenario for perturbative unification is viable; in the sense 
that it is not compatible with the experimental data on a 3 (m^) and the naturalness bound 
on the susy spectrum. 

The question still remains of the corrections due to the heavy degrees of freedom associ- 
ated to the grand unification group (in our case G = SU(5)). As we have seen previously, 
light thresholds introduce non negligible differences in the evolution of the running couplings, 
and the same will take place with the heavy masses (Mj), mainly when the ratio fi/Mj begin 
to approach to one. The evolution of the effective couplings, including only the light degrees 
of freedom, is valid up to scales fi <C Mj, where we are sure about the decoupling of the 
heavy masses Mj. The relationship between the couplings a^ip) and ctQ 1 ^) is obtained 
integrating out the heavy degrees of freedom from the complete action S[G], and is given 
by: 

a^{fi) = aa l (n) + M/x, Mj) . (4.8) 

The functions \{fjL,Mj) include logarithmic contributions of the heavy masses, as well as 
a constant term due to the light degrees [|]]. As we point out in Sect. 3, the logarithmic 
contribution is what we get when aplying the limit Mj 3> // to the threshold functions. 

For non-susy SU(5), we have three basic parameters to describe the heavy masses: My 
(gauge bosons), Mj, (colored Higgs) and Ms (heavy scalar in the adjoint representation). 
For susy SU (5) we have to add the associated susy partners. If we assume, as it is common, 
that the members of each heavy supermultiplet are degenerated, we get0 

, in 66 96, My 4, M$ 20 , M s , A n . 

(4.)A l(ft Mj) = - + -In— - - In— - - In — , (4.9) 

41 Mv My 

4tt Aafo Mj) = — + 8 In — - - 8 In — ^ , (4. 10) 

6 n My 

These expressions are valid at one-loop order for scales \x such that <^ \i <^ Mj. Since 



9 The expressions given in Ref. |2l| for \^{p,^Mj) and also In My included minor mistakes that 
we have corrected here (see Eq. ( 4.1 1| ) and ( |4 . 1 3|) ) . We have also corrected the constant term due 
to the light degrees of freedom that appears in Aj. In Ref. [^9| we reproduced the terms given by 
N.-P. Chang et. al. in EJ. Here, we have recalculated them to be compatible with our definition 
of effective couplings. Anyway, the numerical differences are small, and they do not affect the 
conclusions. 
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heavy masses are typically of order 10 16 GeV, and light masses are expected to be less than 
1 TeV, we choose \i = 10 7 GeV. 

As before, from ([4.8[) we get a system of equations, where now the unknowns will be 



-i 



a G ~, a 3 *(m z ), and In My instead of \nM X - We get for a 3 1 (m^) and In My the following 
expressions: 

^{mz) = -x (3a 2 _1 (m z ) - a^ l (m z )) + ^- Y^F^m;, //) 



3 M$ 3 
In , 

57r fj, 5n 

3tt 



(4.12) 



In M v = — (a x x (m z ) -a 2 1 (m z )) + — ^^F(m;, fi) 



3 , M* 1 , M s 191 
H In — in 



40 // 



// 320 



(4.13) 



where, as before, the sums run over all the light masses. These expressions are valid for 
both effective couplings and step-function approximation. In the first case, the Fj(mj,/i) 
are the approximated threshold functions defined in (|A33|) ; in the second case, we simply 
take Fj(mj,/i) = ln(m 2 //z 2 ) when rrii > mz, and Fi(rrii, fi) = ln(m|/yU 2 ) when rrii < rriz- 
Notice that for a^ 1 \ e and In My \ 9 , the dependence on the scale \i is exactly canceled out 
at one-loop order. For the effective couplings, this dependence remains negligible if we 
maintain the condition rrii "C \i <C Mj. 

The qualitative behavior of a^" 1 (m^) and In My with susy masses is the same observed 
without taking into account the heavy thresholds. 

Our interest now shifts to the heavy mass parameters M$ and Mg. 

Since the value of a^ l {mz) only depends on M$, the limits on this parameter will put 
bounds on a^ l {mz)- The lower bound on Mj> comes from the experimental limits on proton 
decay via dimension five operators. The lifetime for the dominant mode is given by |f27f|: 



r(p -> K + u^) = 6.9 x 10 31 



0.003 sin2p H M$ 



10" 



(3 l + y tK 10 17 f(rriq, rriq, ma) + f(m 9 , m h m$) 



yrs, (4.14) 



where we have introduced three more unknown parameters: the hadron matrix element 
parameter (3, which ranges from 0.003 to 0.03 GeV 3 ; the ratio of the vacuum expectation 
values of two Higgs doublets, tan/3#; and the parameter y tK , which represents the ratio of 
the contribution of the third generation relative to the second generation to proton decay. 
To allow an M$ as low as possible, we take j3 = 0.003 GeV 3 , sin 2/3h = 1 and | 1 + y tK \= 1. 
The experimental limit for this mode is r(p — > K + i? fl ) > 1.0 x 10 32 yr |2~1 
into the following lower bound on M$: 



and this translates 



M $ > 1.2 x 10 20 (f(m 



)) + /(m ? ,m,-, rriw)) 



M«, min 



(4.15) 



The functions /(m 1 ,m 2 ,m 1 j)) come from the dressing of the dimension-five operator with 
the wino exchange, needed to convert it into suitable four-fermion operators for proton 
decay pi ]. With the parametrization we have adopted for susy masses, the combination 



that appear in (|4.15| ) is given by: 
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f(m§, nig, m a ) + f(rriq, m b m a ) = 

1 ^° + 13 ' 5 ln(eo + 7) - ^|ln(eo + 0.5)) . (4.16) 



6.5mi /2 V £o + 6 6-0.5 

So that, M$ mm depends on mi/ 2 and £o, decreasing with both of them. 

On the other hand, the upper limit on M$ and Ms is derived by requiring that the 
Yukawa couplings involving these fields remain as perturbative couplings below the Planck 
scale. This leads to the conditions Mj, < 2My and Ms — 1.8My ||27|1 , which combined with 
(EH) give: 

lnM$ < (a^K) - a^imz)) + Y^g'i{m h /J,) 

5 191 40 

In M s + — ln2 = lnM™ ax - (4.17) 

37 80 37 * V ; 

The upper limit on M$ depends not only on the susy masses but also on Ms. 

We have already all the ingredients (Eq. (|4.12 ), (|4.15|) and (^4.17|) ) to check if the scenario 



of perturbative unification is compatible with all the constraints on a^ 1 (mz), M$ and the 
susy masses. In order to compare, we examine first the results obtained with the step- 
function approximation. In this case, there is no problem in having a^(mz) within its 
experimental range and the susy masses below 1 TeV, as it can be seen in Table |HI[ 

Some comments about these data. The bounds on the heavy and light mass parameters 
are derived from imposing the condition M$ mm (m 1 / 2 , Co) < M^ max {mi, M s ). In this way, 
once we calculated the lower bound on M$ from ( |4.15| ) with the maximum allowed value 
for £ (determined by m™ ax and = 45 GeV), this lower bound on M$ gives us the 

maximum allowed value for Ms. At the same time, the minimum value of Ms gives us, 
on the one hand, the lower bound on £ (upper values on mi/2 and mo), and also gives us 
the upper bound on M$ (Eq. (|4.17| )) and the lower bound on ijiq. These latter bounds 



are derived taking mw2 = 45 GeV and m^=m^, and therefore are independent of any other 
constraints on the susy masses. The same does not occur with the other bounds, which 
depend on the value of the upper limit for the susy masses being considered (mainly m™" 1 ). 
(We have also taken m t = 91 GeV and rrih = 60 GeV to get the values of the Table, except 
to derive a 3 (m z )^ mm ' ) ). 

In principle, we do not have any constraint on Ms mm , except the requirement that there 
is no large splitting between heavy masses. As M$ and My will be around 10 16 GeV, we 
have taken M s mm = 10 13 GeV to give the results quoted in Table |TTT|. With this choice, 
we see that when imposing the naturalness bound on the susy masses, we get mu2 < rriz, 
near its lower experimental limit, and also narrow ranges for m , M$ and Ms: m w 1 TeV, 
M$ w 10 16 ' 7 GeV and Ms ~ 10 13 GeV. We can increase the ranges of rai/2 and m~ h if we take 
a lower value of M E . For example, with M s mm = 10 10 GW we obtain: 45 GeV < mi/2 < 
610 GeV, 415 GeV < mo < ITeV. In any case, if we maintain the naturalness bound of 
1 TeV for the susy masses, it is not possible to maintain all heavy masses of the same order. 
Since the value of M$ is practically fixed by the experimental limits on proton decay, and 
My diminishes with Ms, the most favorable situation will be: Ms < My < Mj>. 

When we calculate with effective couplings, the values of a^ l {mz) are lower than those 
obtained with the step-function approximation, mainly due to the thresholds of the massive 
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gauge bosons. An now, to satisfy the constraints on a^ijnz) and Mq> we need to have m , 
m~ h or mil beyond 1 TeV . The results derived in this case are presented in Table [IV]- We 
see that the lower £ (i.e., the lower squark and slepton masses), the higher will be M$, and 
therefore higher values of m~ h or ran will be needed to get a^{mz) < 0.125. If squarks and 
sleptons are below 1 TeV, the Higgs and higgsino masses have to be beyond this bound. 
Notice that this conclusion is not affected by the value adopted for M^ mm . It is the value 
of M$ mm and the condition a^{mz) = 0.125 which fixes the lower bound for m~ h = ran, 
independently of M s . 

As before, the bounds on the susy masses depend very midly on the values of m t and 
rrih, decreasing as m t and rrih increase. Thus, we have taken the values m t = rrih = 200 GeV 



to give the results in Table |V 



We have seen that, in the context of perturbative gauge coupling unification at one loop 
order, the use of effective couplings and the requirement of having susy masses not too high, 
favor a value of a 3 (m z ) near its upper experimental bound, and this is independent of 
whether or not we include the heavy threshold effects. However, when we try to incorporate 
corrections at two-loop order, or when we take into account heavy threshold effects, the 
maximum allowed value for a 3 {mz) is not enough to get unification, and we have to increase 
the susy masses beyond their naturalness bound. Therefore, when we use effective couplings 
is not so easy to maintain the naturalness bound for the susy masses, as it was when we 
treat the thresholds with the step-function approximation. 



V. COMPARISON WITH OTHER APPROACHES AND CONCLUSIONS. 

In this paper we have studied how threshold effects associated to massive degrees of 
freedom modify the evolution of gauge coupling constants with scale; we have done it for both 
the Standard Model and its Minimal Supersymmetric extension. As a explicit application, we 
have considered the impact of these effects on the perturbative unification of the couplings. 

To begin with, we have examined and specified renormalization procedures for computing 
the threshold functions. In dealing with gauge coupling renormalization, we can distinguish 
two classes of contributions: those coming from gauge boson vacuum polarization, which are 
well and uniquely defined, and contributions coming from vertex and external legs, which 
not so clearly defined. In order to get the latter contribution, we have required that the 
procedure respects the Slavnov-Taylor identities. This will guarantee the universality of 
the renormalized coupling and, at the same time, will constrain to some extent the class of 
particle masses that can contribute to the vertex function T: only gauge boson masses. All 
other massive matter fields contribute to the effective couplings only through the transverse 
component of the vacuum polarization. 

For the QCD gauge coupling, since gluons are massless, the vertex function is simply 
given by the standard term calculated using an MISP. For the SU{2) L x U(l)y gauge 
couplings, the symmetry is broken at low energies, and we have massive gauge bosons. To 
get the explicit expressions for the vertex function, we follow the arguments given in Ref. 



17| . In this way, at the same time we get universal and mass dependent couplings, so that 
in a single stroke we solve another problem related to the electroweak symmetry breaking: 
the misdiagonalization of the neutral mass matrix. 
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We remark that we follow Ref. [17j only in order to define the vertex function T. Our 



effective couplings differ from their definitions. They define the effective couplings including 
only vector boson self-energies, that contain at least, one exactly-conserved matter current 
from the unbroken subgroup SU(3) C x U(l) em . This condition imposes severe constraints on 



the kind of matter fields appearing in the effective couplings there are no contribution 
from neutral and colorless matter fields. Therefore, any relation of the couplings relevant 
for gauge coupling unification will be independent, for example, of the neutral higgses and 
higgsinos. When we study unification with this choice of the effective couplings, we are 
interpolating between the unbroken SU(3) C x U(l) em and the unification group G, but 
forgetting about the electroweak symmetry breaking in the matter sector of the model. 

One can understand the definition of the g'(q 2 ) and g(q 2 ) given in Ref. as a direct 
splitting of the electromagnetic coupling e(q 2 ), such as: 



e 2 (q 2 ) e 



+ ^)+( n, ^(9 2 ) + n^(g 2 )) 



g'\q 2 ) + g\q>) ' (5 ' 1} 

where the n^'s include the transverse vacuum polarization and the appropriate vertex func- 
tion. We choose, instead, to define independently the couplings e 2 (q 2 ) and g 2 (q 2 )- Thus, the 
relation (|5.1| ) only fixes the coupling g' 2 (q 2 ). With our choice, all the particles that interact 
under SU(2)l are present, and affect the evolution of the effective couplings. Of course, the 
high energy limit of the effective coupling is the same for both definitions, as it is dictated 
by the RGE and the decoupling theorem. 



In some sense, our approach is more closely related to the work of Ref. [[12] on light thresh- 
old effects. These authors related the corrections (Sgf) to the low energy gauge couplings due 
to light thresholds (new physics different from the Standard Model) to experimental quan- 
tities measured in scattering processes. Thus, they define independently Sg 2 and 5(g' 2 + g 2 ), 
the first being related to the transverse vacuum polarization of the W^, as we have done. 

Our work agree with the one presented in Ref. |12|] on the importance of light and 



heavy thresholds for the predictions derived from perturbative unification. In spite of this, 
the philosophy of our approach is in some sense different. Taking as unifying group the 
supersymmetric extension of SU(5), they extrapolate the relation given at the unification 
scale, down to the low energy scale mz as 

I2a 2 l - 5«r X - 7a^ \ mz = ; (5.2) 

We only get this relation when using a mass independent renormalization procedure. The 
experimental values of the couplings are extracted on the bases of the validity of the Standard 
Model, but this relationship is obtained when we assume an extended model. Therefore, if 
one wants to make use of the experimental data on the couplings, it is necessary to include 
the corrections due to the new degrees of freedom^. This is also the main issue of some 



10 We would like to rectify the comment on the results of Ref. [12] that we included in the Note 
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recent work [|T^]. In these papers, a calculation of complete susy contributions is carried 
out in order to determine the value of sin 2 Qw{ m z) \jjs f rom experimental quantities, but in 
the context of the MSSM. These threshold contributions evaluated at the scale mz tend to 
decrease the value of sin 2 ^(m z ), and as a consequence there is an increase in the predicted 
value of a^{mz) from unification, hardly compatible with the experimental one. In addition 
to this, we were concerned on how thresholds effects modified the evolution with scale of 
gauge couplings, to understand how in going from Mx to mz these effects will affect the 
relation ( |5.2| ). We have considered a more accurate description of these effects, beyond the 
leading logarithmic correction given by the step-function approximation. When including 
the complete threshold function, not only susy thresholds have to be considered, but also any 
threshold effects due to the massive standard particles. For example, the standard-gauge 
boson threshold effects, which become relevant in realizing or not a perturbative scenario 
for the unification of the couplings, as we have shown. 

We have found that the predictions derived by imposing perturbative unification of gauge 
couplings in the MSSM depend on the procedure chosen to treat the thresholds. For example, 
the effective couplings always favor higher values of a^ 1 (mz) and susy masses, than the 
calculation with the step-function. Moreover, although both procedures give good results 
at one-loop order without including heavy thresholds, when we improve the calculation (2- 
loop order or include heavy thresholds) the unification with effective couplings is in conflict 
with at least one of the constraints we impose on the model: the experimental values for 
ct^niz) and the proton life-time, or the theoretical naturalness bound for susy masses. In 
order to maintain the experimental constraints, we would need susy masses beyond the 
upper bound of 1 TeV usually required because of naturalness reasons. 

At least, this is the general conclusion within the framework of SU(5) unification, where 
heavy threshold corrections to the predicted value of a 3 (m z ) are positive. One way to avoid 
this result would be to considere a different unification group, which gives the reverse sign 
for this contribution [IT3J. We would have to look for a cancellation between "light" and 



"heavy" threshold effects. 



added in proof of Ref. [|29[|. There, we misunderstood the notion of "vertex correction" used in Ref. 
[12]. Of course, the vertex correction due to new physics (susy masses) are not universal, and do 
not have to be included, as they correctly did not [p5| . 
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APPENDIX A: 

In this Appendix we give the analytical expressions for the functions Il^(g 2 ) and T(g 2 ) 
introduced in Section f|. These expressions have been calculated using dimensional regular- 
ization, and can be written as combinations of the following integrals: 

B (ai,a 2 ) = / dx\n(aix + a 2 (l — x) + x(l — x)) , (Al) 
Jo 

B 3 (ai,a 2 ) = / dxxil — x) \n(aix + a 2 (l — x) + x(l — x)) , (A2) 
Jo 

rl rx 

C (ai, a 2 , £a 2 ) = dx dy\n(ai(l — x) + a 2 (x — y + £,y) + x(l — x)) , (A3) 
Jo Jo 

A(Gl ' G2 ' ea2) = Io dX Io dy 'Mi -x)+ a*{x-y + £y) + xjl - xj) ' (A4) 

Ds(ai ' G2 ' eG2) = /o^ /o^ fa (1 - x) + Jx-y + ft/) + sfl - s)) ' (A5) 



;i-0^(ai,a 2 ) 



(ai(l - x) + a 2 (x - y + £y) + x(l - x)) 
1 /"j /"*, , ^{l - x + i{x - y)) + a 2 y + y{l - y) 



H** A/ 1 " T = x ~ y T y +y yr y v (A6) 

2a 2 ,/o 7o ai(l - x + +£a 2 y + 2/(l -j/) 

and the linear combinations: 

D+ =D l + D, , Dii = D l - D 3 . (A7) 

2 

Here, we have introduced the variable a« = — l —, where (— q 2 ) is the euclidean momentum. 

—q 2 

The general transverse contributions to the 11^, which depend on the kind of masses 
running in the loops, are given by (£ is the gauge parameter): 

(a) Gauge Boson + Ghost: 

( 4^)( 0l , 02 )=-i{(f-,)g-l„^) + ? 

-2S (ai, a 2 ) - 10JB 3 (ai, a 2 ) + 2B 3 (£a 1 , £a 2 ) 
-(1 - (l + C (a u a 2 , £a 2 ) + C (a 2 , a u fai) 

+ (1 + ai)D^(ai, a 2 , £a 2 ) + (1 + a 2 )Df 3 (a 2 , a 2 , £a 2 ) 
+(l-£) 2 #(a 1 ,a 2 )j. (A8) 

(b) Scalar-Gauge Boson: 

(4 7 r) 2 n( SG )(a 1 ,a 2 ) = -(1-Oa2^(ai,a 2 ,eo2), (A9) 

Here m\ is the scalar mass, and m 2 the gauge boson mass. 

(c) Scalars: 
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(4n) 2 rt s \a u a 2 ) = ± |0- In -J^ -W (a u a 2 ) + 12 B 3 (a u a 2 )^ . (A10) 
(d) Fermions: 

(47r) 2 n^(a 1 ,a 2 ) = ^{f--ln^ -6fi 3 (ai,a 2 )| . (All) 



3 [ y e fi 

In the expressions above, we have not included the couplings or group factors. The 

scale "it" is the unit of mass introduced in dimensional regularization, and e is given, as it is 

2 2 

usual in Modified Minimal Subtraction, by - = 7 + In Air. 

e n — 4 

In each particular case, we get: 

nL(? 2 ) = e 2 \2rt G \a w , a w ) + U^ SG \^a w , a w ) + U^(^a w , £a w ) 

+ £G}n<'Ha/,a / )}, (A12) 



V-wwiv 2 ) = g 2 {2s 2 w rt G \0, a w ) + 2c 2 w rt G \a z , a w ) 
+rt SG \a h , a w ) + s 2 w rt SG \£a w , 0) + s 4 w ^ SG \^a w , a z ) 

+ ±rt s \a h ,Za w ) + ±rt s \Za z ,Za w ) + ± £ n^(a /l; a /a )] , (A13) 

W) = ^ 2 {3n^(o,o) + \ £ nW( B/)fl/ )} (am) 

quarks 

In order to obtain the function T(g 2 ) = — (g 2 + g /2 )n§ A (g 2 )/(gg')m|, we need the lon- 
gitudinal term of the diagrams in Fig. (|9|). The diagrams (^.a) and (||b) are related and 
we do not need to calculate both explicitly; the same kind of diagrams contribute to the 
longitudinal vacuum polarization of the photon, which we know to be zero. Therefore, if 
we name A(a w ) the longitudinal term of ([|a), and B(a w ) the corresponding to (§.b), they 
must verify the following identity, 

^aa(Q 2 ) = e 2 m 2 w (A(a w ) + B(a w )) = -> A(a w ) = -B{a w ) . (A15) 

So that, for Yl^ A (q 2 ) we obtain: 

g 2 g' 2 

^■za(q 2 ) = -e , =m 2 v A(a w ) + eg— — -m z m w B(a w ) = gg'm 2 w B(a w ) , (A16) 

V9 +9 g +g 

where, 

B(a w ) = {^y^ ^ - In ~ 2B (£a w , a w ) 

+(1 -QC (£a w ,a w ,£a w ) - 2(1 - £)D^(£a w , a w , £a w )\ , (A17) 
And finally, we get for T(q 2 ): 
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771 

T(g 2 ) = (g 2 + g' 2 )^fB(a w ) = -g 2 B(a w ) . (A18) 
As we did in the main body of the paper, on the following we redefine H"^ = gfUfi and 

I' /rr. 

Now, we write the expressions for Uf^q 2 ) + 2Ti(q 2 ) in the Landau gauge, and separate 
in each term the divergent part: 

(An) 2 (Tl T AA + 2T(q 2 )) = 

13 (2 -q 2 \ (2 -q 2 \ 

— - ( - - In — - + F G (a w , a w ) \ - 3 ( - - In — — + F r (a w ) J - F SG (0, a w ) 

+\ - ^ ±£ + ^(°» °)) + \ E Q} - ln ^ + ^ °/)) ' ( A19 ) 

(47r) 2 (n^ + 2r(g 2 )) = 

13 /2 -q 2 \ (2 -q 2 \ 

— g- ( - - hi + s^F G (0, a w ) + c 2 w F G (a z , a w ) \ - 3 f - - In — + F r (a w ) J 

1 (2 -q 2 \ 

-F SG (a h , a w ) - s 2 w F SG (0, 0) - s^F SG (a z , 0) + - 51 ~ ~ ln — + ^V(a/n a / 2 ) 

d doublets \ £ A* / 

+^ - In ^ + ^ fl (a fc , 0) + ^F s (0, 0)) , (A20) 
(47r) 2 (nj s + 2r 3 (g 2 )) = 

-ll[H_i n ^ + C 7te]+| £ ^-ln^ + F^a,)) . (A21) 

\ " / quarks \ r 1 / 

The functions Fj(a,i) contain the threshold effects associated to massive degrees of freedom. 
Each of these functions has the property that they tend to a constant in the limit where the 
masses can be neglected as compared to the scale, and they behave as lna^ + 0(1/ di) in the 
limit of heavy masses compared with the momentum scale; except the function F SG (ai,a,j), 
which tends to zero in both limits. Having in mind these limits, we can study the behavior 
of the effective couplings, 



v 2 =ql 



for several ranges of scale: 
(a) Limit q 2 > ql > m 2 : 



(An) 2 


(An) 2 


e 2 (q 2 ) 


e 2 (ql) 


(An) 2 


(An) 2 


g 2 {(f) 


g 2 (<i) 


(An) 2 


(An) 2 


gi(Q 2 ) 


g 2 M) 



13 4 
1 ~3 + 3 



V 3 6 3 doMetsJ Q 2 9 2 {%) 6 q 2 

+ ? EV4^-71n4- (A24) 
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Thus, we recover in this limit (negligible masses) the usual expressions derived by minimal 
subtraction. 

(b) Limit q 2 < q 2 <C m|: in this case, the massive degrees of freedom must decouple 
and not contribute to the effective couplings. When working in the Landau gauge, we have 
to pay special attention to the contribution of the Goldstone bosons and ghosts. In the 
Landau gauge, we treat these degrees of freedom as massless, and we could argue that the 
limit of heavy gauge masses will not affect them. If we proceed in this way, we will get 
the non-decoupling effects proportional to the gauge masses. The correct result is obtained 
reverting the limits: first, we take the limit m 2 /q 2 — > oo in a general gauge, and then we 
take £ — > 0. We will assume this procedure whenever we take the limit of heavy masses^]. 

(c) Limit q 2 , <3Zi m 2 <^ <f"- m this case we recover, in part, the results obtained by using 
the step-function approximation. For example, for the electromagnetic coupling we get: 



^ = ^+ (llJ(q 2 ) + 2T l (q 2 ))- lim (uj(q 2 ) + 21^ 

- (4?r)2 -l(\n m2w 




e 2 (q 2 ) 

In the last line, we simply have absorbed the constant terms in the logarithms, except for 
the factor "2/3", which has a different origin^. The differences between this expression and 
the one obtained by the step-function method are in these constants. Their origin is better 
understood if we study the (3 function directly, which is given by: 

= -j/ 9 (%) - 3/rM + ~/ s (0, 0) + ^ E Q}ffM > (A26) 

where we have defined, 

13 , , , 13 / dF G ( )\ dF SG {0 

-jf 9 {a w ) = -\l d ^ q2 j - , (A27) 

/rW = 1 "7V' (A28) 



11 For example, we have considered this when we talk before about the limits of the function 

F SG (ai,aj). 
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It is a typical threshold effect of the massive gauge bosons, produced by the regularization 
method employed. Indeed, it is not present when using dimensional reduction 36]. 
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We can seen in Fig. (|10D , where we have plotted //(a/), that the behavior of the functions 
fi{dj) is similar to the step-function. When we integrate (3 e from q 2 , to q 2 , we can divide 
the integration range into [q 2 ,m 2 ] and [m 2 ,q 2 ], and approximate the functions /j(aj) in 
each interval for the corresponding limit [m\jq^ — > oo, m 2 /q 2 — * 0). In this way, we 
reproduce exactly the results of the step-function approximation. Nevertheless, if instead 
we use the Taylor expansion of the functions fi(aj) to first order in Oj (improving the latter 
approximation), we would get for a fermion of mass m/, 

ri 2 j" m2 /i 2 1 rm 2 f/q 2 

- /, ff(m 2 f /y)dliay ~ - / __<nny + / (1 + c%)d\ny 

9 -i 2 2 2 

-^2-^ + ^-/ + 4-/-cJ = lnc-l, (A30) 



cymj Cj q" J q* J q" 

where once again a constant appears in the logarithm. But, neither of the two approxima- 
tions, step-function or Taylor expansion, is valid, because when we integrate we cross the 
region m 2 /q 2 ~ 1, where none of these approximations is defined. Due to this, the best 
procedure is to approximate the complete function, //(a/), by an expression valid for the 



full energy range, as it is |37] 



m^tt^t,- (A31) 



c f a f 

When we integrate with this function, we get: 

- f * f f {m)/y)dlny * - [*"''* -J— dlny = Inj + C/ j . (A32) 

Therefore, we have seen that, for any massive particle, its contribution to the effective 
coupling can be approximated by a logarithmic function of the form: 

L k (q 2 ,m 2 )=ln q l + Ckm2 2 . (A33) 
q z + c&m z 

With this function, we recover easily the limits (a) and (b) for the effective couplings and, on 
the other hand, we get a more precise behavior of the thresholds for intermediate scales. The 
value of the constant depends on the kind of particles we treat (fermion, scalar,..), and if 
it is or not degenerated with its partner in the loop. In particular for the Standard Model, 
with or without supersymmetry, and the range of energies of interest, the contributions from 
loops of non-degenerate particles can be treated as completely degenerate, or that one of 
the masses can be neglected with respect to the other. A numerical study shows that the 
best values of the constant for these cases, are those given in Table |T[ 

In the following, and to close this Appendix, we include the analytical expressions of the 
functions Fj (01,02): 

5 1 

F F (ax, a 2 ) = - - - lnaia 2 - 2(a t + a 2 ) - 2(a x - a 2 ) 2 

--(ax - a 2 )(3(ai + a 2 ) + 2(a x - a 2 ) 2 ) In — 
2 a 2 
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+-(1 - oi - a 2 - 2(ai - a 2 ) 2 )^i2^i2 , 
8 1 

F s (ax, a 2 ) = - - - In aia 2 + 4(ai + a 2 ) + 4(a x - a 2 ) 2 

+-(ai - a 2 )(3 + 6(ai + a 2 ) + 4(ai - a 2 ) 2 ) In — 
2 a 2 

+-(1 + 2(ai + a 2 ) + 4(ai - a 2 f)R 12 L 12 , 

^sg(«i, a 2 ) = ~(a 2 - 2a 2 + 4a 2 ai 
o 



«2 



+ (1 + 3ai + 3ai(ai — a 2 ) + (ai — a 2 ) ) In — 



a i 



1 + cii, 



1 + 2ai - a 2 + (ai - d 2 ) 2 )R l2 L 12 - 2(1 + ai) d In ' — ) 

ai 

2 1 + a 



F r (a) = --lna+^- — ln(l + a) - (1 + a + ^-) In 

o o da 3 a 

3 

— F G (ai, a 2 ) = A(ai, a 2 ) + ai_B(ai, a 2 ) + a 2 £?(a 2 , ai) 

C(ai,a 2 ) C(a 2 ,ax) a 2 . ax E(a 1 ,a 2 ) 
H 1 1 ^(ai, a 2 ) H -L>(a 2 , d\) H 



(A34) 



(A35) 



(A36) 
(A37) 

(A38) 



A(ai, a 2 ) 

B(ai, a 2 ) 
C(oi, a 2 ) 
D(a 1; a 2 ) 

E(ai, a 2 ) 



121 



7 13 a 2 

- In aia 2 + 6aia 2 + — aia 2 (a 2 — ai) m — 
3 3 ai 
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+(— + 3aia 2 )i?i 2 Li 2 , 

„ 8 .3 17 7 2 ., ai .13 4 , 

-4 - -d x + - - —ax - -aj In — + -a^R^L 

3 8 o o a 2 12 3 

7 7 1 

-— In d X d 2 + — ln(l + a 2 ) + --Ri 2 L i2 , 



12 , 



1 

24 



(23 + 17a 2 - 2a 2 - Ad\) (2 ln(l + a 2 ) - In d X d 2 ) 



+ (15 + 2a 2 - Ad 2 2 )R 12 L l2 
1 
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(In did 2 - 2 ln(l + oi) - 2 ln(l + a 2 ) - -Ri 2 L 



12/ 



where we have introduced: 

-R12 = 
L 12 = In 



1 + 2(di + d 2 ) + (ai 
1 + di + d 2 — R 12 



«2 



1 + di + d 2 + i? 
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(A39) 
(A40) 
(A41) 

(A42) 
(A43) 



(A44) 



We note that as regularization procedure we use dimensional regularization to work with 
the Standard Model, and dimensional reduction to work within the MSSM. Contributions 
from fermion and scalars are the same in both procedures, the only difference being in 
the gauge boson contribution. For translating to dimensional reduction, we have only to 



subtract a constant term "2/3" from equation (|A39|) . 
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APPENDIX B: 



In this appendix we approximate the effective couplings at the 2-loop order. The RGE 
for gauge couplings at the 2-loop order are given, in general, by ||38|,|7| 

When we calculate the coefficient bi and 6y using an MISP, we obtain constant values for 
bi and 6^; but when using an MDSP and we include the dependence on the masses, these 
coefficients gain a dependence with the scale /x, through the ratios (m^///). The take 
into account the threshold effects at 1-loop order, that we have calculated in Appendix 
A. Now, we want evaluate which is the threshold contribution at 2-loop order, included 
in the with respect to that at 1-loop order, without actually performed the explicit 

calculation. 

When we integrate ( |B1| ) from mz to an arbitrary scale /i, we get: 

1 /"A* bi(uf) 

or 0) = a t ~ Qm) \ x _ loop -— * dlna^n') . (B2) 

In order to evaluate the contribution of the second term, we choose an intermediate scale, 
m z < H\ < /i, such that <C fi±, for all the masses m ; . For example, for the Standard 
Model and its supersymmetric extension, we have < 1 TeV, and thus it is sufficient if we 
take n\ = lOTeV. In the range \i! > /ii, the thresholds at 1-loop order and 2-loop order 
are negligible, and then the coefficients £>»(//') and remain as constants. Therefore, we 

can write: 

1 rWK lT _ ,. A i rk 



-i/ % i / \ i 1 f^ 1 bij{u ) , 1 f 11 on , 

Ma 



--r'Mi— + 4^'"w + ? Ay - (B3) 

Now, the 2-loop threshold effects are included in the function Ay, and this is the function in 
which we are interested. To evaluate the order of magnitude of this function, we approximate 
the integrand in A^- by a straight line in the variable In a 3 - (//'), 



'bij(fi') b 



'j 



+ CijlnaAji') , (B4) 



with the condition that it passes through the points mz and /ii, i.e., 

fj! = m z ^ an + c {j In djK) = f ^|^y - ^ = Bij{m z ) , (B5) 
//' = //!-> a ti + In aj(/ii) = ^^^y ~ ^1 - • ( B6 ) 
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With these conditions, the expression for A i3 - is reduced to: 

87T 

The final step consists of taking at 2-loop order the same kind of threshold functions, 
fi{mi/ji), that at 1-loop order, i.e., a function bounded between and 1. For example, 
623 ~ 8 + 4/ t op, and then the value of A23 will be between 0.013 (/ iop = 0) and 0.003 
(/top = 1)- Checking for different values, the corrections to a" 1 , given by Aj = J2j A^, are 
never higher than 0.3%, and therefore, negligible. 

Thus, the 2-loop effective couplings can be approximated by: 

where we only need to include explicitly the threshold functions at 1-loop order. 

This approximation is independent of the chosen regularization procedure, in particular 
if dimensional regularization or dimensional reduction. The regularization procedure does 
not affect the coefficients of the RGE's for the gauge coupling when we work with MS and 
always preserve its form (at one or two loops). With an MDSP, it easily seen that the function 
bi(n) at one -loop order does not depend on the regulator. When working at two-loop order, 
we would have to check how the explicit expressions of the bij (/^-coefficients depend on the 
regularization method. However, here we have only made use of the general properties of 
these functions, i.e., their behavior in the limiting cases mj\i — * and \ijm — > 0. And 
this behavior is independent of the regularization procedure, as well of the particular mass 
dependent subtraction procedure one uses. 
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FIGURES 



FIG. 1. The 1-loop particle contributions to the IlJ g function, (a) Gluon contributions plus 
ghost (uJi) contributions, (b) Fermion contributions. 

FIG. 2. Evolution of the three couplings of the SM calculated with the three 1-loop procedures: 
effective couplings (solid lines), step-function (dashed lines), and MS (dotted lines). 

FIG. 3. Derivative of the effective coupling a^ 1 ^) respect to In fi (solid line). The dashed 
line shows the derivative within the step-function approximation. 



FIG. 4. Evolution of the three couplings of the SM at 2-loop order with MS (dashed lines), and 
1-loop order with a mass dependent method (solid lines). We also plot 1-loop MS for comparison 
(dotted lines). 

FIG. 5. Evolution of the three couplings of the MSSM at 1-loop order calculated with 
step-function (solid lines) and effective couplings (dashed lines), for different values of the susy 
mass parameters: (a) m^ 2 = 45 GeV and tuq = mz, (b) m^ 2 = m$ = 1 TeV; we take 
m~ h = m + = niH = mo and m t = m^ = 200 GeV . We include in the plot the experimental 
error band for each coupling. 

FIG. 6. Same as Fig. 5. a, but now without including the threshold effects due to the massive 
gauge bosons in the evolution of 1 ■ 

FIG. 7. Values of a^ 1 (mz), compatible with the unification condition (Eq. 19), calculated 
with effective couplings at 1-loop order, for different values of = mn- mz, 1 TeV, lOTeV, 
100 TeV. Dotted lines are the experimental limits on a^ 1 (mz) and w.1/2; solid lines are for 
mt = m^ = 200 GeV, and dashed lines for mt = 91 GeV and m^ = 60 GeV. The straight lines 
(solid for m t = m^ = 200 GeV and dashed for m t = 91 GeV, m^ = 60 GeV) are the upper limit 
obtained for a^ 1 (mz) when imposing Mx = 10 16 GeV {m~ h = mn increase along these lines from 
bottom to top). The allowed region for a^ 1 (mz) are to the left of the straight lines, and between 
the dotted lines, 8 < a^ 1 (mz) < 9.2 and > i5GeV. 

FIG. 8. Same as Fig. 7, but with a 3 ' 1 (m^) calculated with the step-function approximation. 

FIG. 9. The 1-loop particle contributions to the I1^ A function, i.e. to T' . (a) Gauge boson 
plus ghost, (b) Scalars and gauge bosons. 

FIG. 10. Function //(a/) respect to —logio(af) = log\o(q 2 /mj). As it is shown in the plot, the 
function tends to 1 in the limit of neglecting mass respect to q 2 (—logio(af) — > oo), while it tends 
to (decoupling) in the limit of heavy mass {—log\o{af) — > — oo). 
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TABLES 



TABLE I. Experimental values of 03 (mz). 



Experiment 


Central Value 


Error 


ALEPH jets 


0.125 


±0.005 


I J Tj 1 j L ILL jeub 


n 1 1 ^ 


znu.uu / 


DELPHI (e+e - ) 


0.118 


±0.005 


L3 jets 


0.125 


±0.009 


OPAL jets 


0.122 


±0.006 


OPAL r 


0.123 


±0.007 


J/* 


0.108 


±0.005 


Deep Inelastic 


0.111 


±0.005 


UA6 


0.112 


±0.009 



TABLE II. Fitted values of c k . 



Massive particles c& 

1 Gauge boson 2 

2 Gauge boson 5 

1 Fermion 2.5 

2 Fermions 5 

1 Scalar 4 

2 Scalars 10 



TABLE III. Values calculated with the step-function approximation^ 



^.max 

m 


Cmax 


M min 


M max 


„,rran 


m l/2 


Cmm 


m mm 


M max 


„,max 


1000 


494 


-LQ16.58 


10 14.5 


0.112 


91 


121 


742 


10 16.79 


0.125 


2000 


1975 


-LQ16.06 


10 16.3 


0.108 


358 


31 


742 


^q16.79 


0.125 



a Mass values in GeV. 



TABLE IV. Values calculated with the effective couplings 21 . 



m max 


cmax 


m min 


M mm 


M max 


~,min 
3 


m max 


tmm 
?0 


M max 


s,,rnax 
3 


1536 


1165 


1536 


-LQ16.26 


^q16.49 


0.125 


45 


1165 


1Q16.26 


0.125 


3000 


4444 


370 


10 15.75 


-^q16.53 


0.120 


239 


158 


-^q16.50 


0.125 


4000 


7901 


198 


^q15.53 


-^q16.54 


0.118 


794 


26 


10 16.61 


0.125 


5000 


12348 


121 


10 15.35 


10 16.55 


0.116 


2784 


3 


10 16.69 


0.125 



a Mass values in GeV. 
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Figure 5. a 




Figure 5.b 
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